CERN-TH/98-309 
FAMNSE-98/37 
LPTHE Orsay-98/60 

hep-ph/98 10449 

Addendum to Finite-size effects on multibody neutrino exchange 

A. Abada 

Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
asmaa.abada@cern.ch 

O. Pene and J. Rodriguez-Quintero 
Laboratoire de Physique Theorique et Hautes Energies*, 
Universite de Paris XI, Bat 211, 91405 Orsay Cedex, France 
pene@qcd.th.u-psud.fr, jquinter@qcd.th.u-psud.fr 

Abstract 

The interaction energy of the neutrons due to massless neutrino exchange 
in a neutron star has recently been proved, using an effective theory, to be 
extremely small and infrared-safe. Our comment here is of conceptual order: 
two approaches to compute the total interaction energy density have recently 
been proposed. Here, we study the connection between these two approaches. 
From CP invariance, we argue that the resulting interaction energy has to be 
even in the parameter b = —Gpn n /\/2, which expresses the static neutrino 
potential created by a neutron medium of density n n . 
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The long-range neutrino-mediated interaction between neutrons in a dense core, such as 
a neutron star, has recently been studied and shown to be extremely small. Following 
the method of Schwinger |7|, the interaction energy W can be computed as 

W = (6\H\6) - (0\H \0) , (1) 

where |0), |0) are the matter and matter-free vacua, respectively, and H (H Q ) is the Hamil- 
tonian in the presence (in the absence) of the star medium. For convenience we will call, 
in this note, "matter vacuum" the neutrinoless ground state in the presence of a star. This 
first step is crucial: the interaction energy is equal to the shift of the zero-point energy of a 
neutrino due to the presence of the star. We have demonstrated PS that there is a non- 
vanishing zero-point energy density difference (W) between the inside and the outside of 
the star, which is due to the refraction index at the stellar boundary and the resulting non- 
penetrating waves. In Ref. ||, this effect was shown analytically and numerically to be the 
dominant one and lead to an infrared-safe total energy density. This result is in contradic- 
tion with the previous claim in || that there must be a lower bound on the neutrino mass to 
ensure the existence of stars. The latter "catastrophic result" is a consequence of summing 
up large infrared terms outside the radius of convergence of the perturbative series. The 
use of an effective Lagrangian which can be exactly solved, allows this non-convergence 
to be circumvented and to provide the result to a good accuracy. The effective Lagrangian 
writes: 

£ e ff = iV L ^ v L {r) - b 77 L 7 z/ L 6(R - r) , (2) 

where R is the radius of the neutron star and b = —Gpn n j \[2 summarizes the static potential 
felt by the neutrinos, which is generated by the neutrons through Z -exchange. The neutrons 
are treated as static, and represented by a uniform axial-charge density (being electrically 
neutral the neutrons interact with the Zq only the via axial current) while n n is the star neu- 
tron density. Notice that the potential is attractive, with depth b, for neutrinos, which then 
condense, and repulsive for antineutrinos. It is important to stress that this effective La- 
grangian, Eq. (0), is non-CP-invariant since it describes the neutrino in a non- CP-symmetric 
background: the neutron star. In fact, it is easy to see that CP C{b) (CP) -1 = C(—b). 

As a consequence of the effective theory used, expression ([!]) is formally an ultraviolet- 
divergent quantity, which needs to be regularized. This can be written as 

W = Y,E l -El (3) 

where i runs over the negative neutrino energy levels. Another method consists in the 
following symmetrization: 
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W sym = -1/2 - E°) - J2(E t - E?)) . (4) 

\i>0 i<0 / 

As a matter of fact, Eqs. (|3p and (|4j) could be considered as alternative definitions of the zero- 
point energy difference, i.e., of the vacuum energy for the theory [0. Both definitions are 
equivalent, as we shall see, only if the effective theory is symmetric under CP transformation. 

Using the latter expression, the result is found to be even in the parameter b ||||. Using 
the former one, we had found it to be odd H and the authors of @|| stressed that this 
last result was incomplete. Their criticism triggered the present study which lead us also, 
although with different arguments, to the conclusion that the result should be even in b. 

To summarize we will mainly discuss two issues: the proper description of the vacuum 
energy and the UV dependence of the result in order to clarify the differences between 
the approaches in 0|§-||, which, needless to repeat, all agree on the main issue: massless 
neutrinos do not imply any catastrophe. 

The energy density w(x), the integral of which is W in Eq. ([]]), is given in terms of 
Hamiltonian densities by 

w(x) = (6\n(x)\6) - (o\n (x)\o) . (5) 

In order to compute it following Schwinger JF], we write it as: 

} ; (6) 

the index (0) refers to the free vacuum and Sf is the usual propagator, 
Sf(x, y) 7o = 9(x° - y°) t lMS)V4U(fl 

- %° - x°) t M^IAV) e~ lE{x °- y0) , 

y~K,E<0 

(7) 

where Hip KyE = Eip Ki E, k summarizing the other quantum numbers, for instance the mo- 
menta. From Eq. (|7|), w(x) can be rewritten using the notation k ,e = J2kI dE in the 
following two ways. 

i) The first one consists in taking the limit y° — > x° with y° > x°: 

^Tk,E<0 

~t E^{x)^ E {x). (8) 

7k,E<0 

Obviously this choice is equivalent to @. 

ii) The second one corresponds to taking the symmetric average of the limits y° — > x° 
with y° > x° and y° — * x° with y° < x° in Eq. (^), as done in [Q and also in B, which gives 



w(x) 



—% 



tr 



\7o 



S F (x,y) - S { p'(x,y) 
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^symV^J 



0^ Ei)l E (x)^ K>E (x) 

l ^Tk,E<0 

1 ^ M ( §(5)A(5) 



,£<0 



4i muw^w 



(9) 



corresponding to (H). The symmetrization over the limit on the time components is equiv- 
alent to the symmetrization over positive and negative energies as a consequence of time 
boundary conditions of the usual Feynman propagator^. Moreover, this symmetrization 
leads to an even behaviour in b: if one takes into account the transformation of the La- 
grangian under CP, CPC{b){CP)~ l = C(—b), Eq. (|9]) can be rewritten as 



w sym (b, x) = - T E Ul tE il) K ,E(b, x) + 
ipl,E^K,E{- b i x) ~ 2^1,^,i?(0, x)} ; 



(10) 



where the dependence on b of the eigenfunctions is explicitly written. 

In Refs. [|]|| this result has been confirmed, while by using the non-symmetric expression 
(|) the energy was linear in b ||. In order to understand this difference, we can subtract 
Eq. © from (g): 



w{x) - w sym {x) 



fi d 

\ 2 dxo Tk,e 

-^E ] m%{z 

i d 



}y-* 



-iE(x -yo) 



2 dx 
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-(o\H {0) (y)MS\x)]\o) 



(ii) 



y— >x 



where [W, stands for the canonical anticommutation relation of the fermionic field: 

Jk,E>0 

+ t. ^,E(^ iEx °dl E (12) 



T7i'E<0 



It is worth pointing out that the use of momentum-space Feynman propagators to rewrite the 
energy density, w(x), implies an implicit symmetrization over positive and negative energies because 
of the Fourier integration over the time and of the time boundary conditions. Thus, in Feynman's 
picture the symmetrization is assumed. 



4 



the "coefficients" b K ^E and d K E being the usual fermionic annihilation and antifermionic cre- 
ation operators; k' is obtained from k by performing the canonical transformation generating 
the positive-energy antiparticle states from the negative-energy ones. 

From the CP transformation: ty(xo,x) — > i^Jo^ixo, —x) ||, it is easy to see that 
the last r.h.s of (JTTJ) changes sign under CP. Therefore ( ]TT| ) vanishes if the vacua ( |0), |0)) 
are CP-invariant. An equivalent way of seeing this is to remark that in a CP-symmetric 
vacuum, there is an (E, k — > —E, /^-symmetry, inducing a cancellation in the sums of the 
first r.h.s in (|TTJ). It results that in a CP-invariant vacuum, ([|) and ([|) or equivalently (Q) 
and (j^) are equal. 

In a non- CP-symmetric vacuum |0), namely a neutronic vacuum leading to the non-CP- 
invariant effective Lagrangian for the neutrinos (Q), which has been used by the authors of 
Refs. Eq. flTT| ) does not vanish. Indeed, the axial charge (to which the parameter b is 

proportional) changes sign under CP transformation |§ and, in fact, our effective Lagrangian 
is only invariant under the product of CP (which changes neutrinos into antineutrinos) and 
the operation that changes b into —b (which changes the neutron star into an antineutron 
st mentioned by Kiers et Tytgat in H). 

Therefore CP invarianceQ imposes the invariance of the "vacuum" energy under the 
exchange b — > —b: the energy due to massless neutrino multibody exchange inside a 
neutron star is the same as the one inside an antineutron star (see Fig. |l|). 

In principle, the problem of matching the effective theory with the underlying one is 
related to the appropriate choice of the ultraviolet regularization. A detailed and faithful 
description of the transition from the underlying theory to the effective theory would lead to 
a UV-regularization, bringing automatically into the effective theory the wanted properties 
such as the CP symmetry. However, this is a very complex task. The limited goal faced in 
Refs. justifies a naive 11 description" of UV physics; for instance, the use of a simple 
cut-off in Ref . 0] . We will hence invoke the general symmetry properties of the underlying 
theory (CP) to constrain the description of the vacuum energy in the effective theory (|j), 
exactly as the chiral symmetry of QCD constrains the effective chiral Lagrangian. 

Eq. flT0D, and hence (^|), is clearly invariant under b —>■ —b transformation. The latter 
equation is no longer equivalent to (||) since the difference given by Eq. (^TJ) is non-zero, the 
Dirac equation obtained from (0) is not symmetric in the exchange E,k — > —E, k,'. Thus, 
the CP symmetry of the underlying theory, i.e. Q.C.D., imposes the choice of (||) in order 



2 In this note we neglect CP violation in the Standard Model. Still, if the neutron stars are 
considered at equilibrium, it suffices to invoke CPT to impose an equal mass to neutron and 
antineutron stars. 
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(a) 



(b) 



n 



n 



(c) (d) 

TT < " IT 

FIG. 1. T/ie diagram (a) is changed to (d) and (b) to (c) by reversing the time arrow. In the 
Feynman's picture, reversing the time arrow implies a change of fermionic lines into antifermionic 
lines. Thus, CP invariance guarantees that the total result will remain the same for both neutron 
and antineutron stars, (a) + (b)=(d) + (c). The dashed double line represents the neutrino propagator 
in the medium of the star. 




to compute the interaction energy density. 

For practical calculations we need to remember that our effective theory is valid only 
under the assumption that the neutrons are static and remain so until the energy scale 
C ~ 100 MeV. Moreover, there are other energy scales that limit our theory, as the one 
related to the confinement of QCD (~ 1 GeV), the one related to the mean free path of 
neutrinos in the star... The physics related to these UV cut-offs has been discussed in 
Ref. |J. The use of a cut-off in the energy integration can be alternatively interpreted as 
the assumption that a repulsive core prevents the neutrons from "piling-up" in space, as 
noticed by Fischbach ||. In our effective theory of Eq. (0), b can be naively interpreted 
as the coherent and homogeneous amplitude for static neutrons to interact with neutrinos. 
The energy cut-off expresses the scale at which neutron recoil, repulsive cores, quark and 
gluon substructures, etc., induce the neutrino interaction with non-static neutrons to become 
incoherent and inhomogeneous. 

In 0, using one of these physical cut-offs, C, we have found for the energies, using 
respectively Eqs. (|) and (||): 

W ~ -bC 3 R 3 . (13) 
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W sjm ~ -b 2 C 2 R* . (14) 

Eq. ([14]) ) which uses (j|), is compatible with the CP symmetry of QCD. Beside the cut- 
off explicit regularization, it is probable that the result should also depend on the way of 
performing the summation over the energy. Indeed, this is the first step in the regularization 
process Q. 

In computing W sym , Kiers and Tytgat M counted the energy levels of the Hamiltonian 
by putting the system in a large box, which was a method different from the one in || , and 
their result was 

W£ ~ -b'R 3 - (15) 

Let us present a few comments to understand the difference between fll4"[) and (|15|). 
Both are even in b and respect CP symmetry. Next, the sums leading to (|i"4"|) and to (|I~5|) 
correspond to different orderings of the same energy levels. Reordering a sum would not 
change the result if we were not speaking of divergent series which have been regularized 
in ways that turn out to be different. This can be illustrated simply by applying both 
summation methods ( || and [|J]) to the (1 + l)-dimensional toy model we presented in 
Ref. 0: we found that W = because there was a one-to-one correspondence between the 
energy levels inside and outside the star; by putting the (1 + 1) star in a box, as done in 
Ref. 0], we verified that the result was not zero, although it was extremely small. So the 
way of summing over the energy levels gives a difference in the final result. 

Furthermore, the main difference between the two symmetrized results, Eqs. ([L4]) and 



(|15|), is the power in the parameter b. The origin of that discrepancy comes from neglecting 
the two-body contribution to the energy density performed by the authors of M . They argue 
that the two-body interacting energy behaves as R 2 , its contribution to the energy density 
scales as 1/R and this becomes negligible in the large i?-physical regime. This seems to 
contradict previous studies that supported the existence of a well-defined neutrino exchange 



two-body potential (10 



- l — . (16) 
4-7T 3 \f\ — r 2 \ 5 

The integration of the latter potential over space is UV- diver gent. Some regularization 
procedures lead to a vanishing result, for example, Pauli-Villars in []l],[|, and dimensional 



3 The classical example is that of the two possible orderings of the non-convergent series (— l) n : 

E„(-l) n = 1 + (-1 + 1) + (-1 + 1) + • • • and En(-l)" = (1 - 1) + (1 — 1) H , which are both 

valid to describe the series. 
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regularization in |4[]. This vanishing of the integrated two-neutrino positive potential ([16]) 
is physically surprising, and is equivalent to assuming a negative distribution located at 
ri = r 2 . Another regularization, such as a simple UV cut-off: \r\ — r 2 | > r c certainly leads 
to a positive non- vanishing result oc \jr\. 

If one still accepts this vanishing |4j] , then no 6 2 -terms remain, and only diagrams with 
at least four neutrons contribute; those contributions start at b 4 . However, in Ref. ||, the 
authors stated that for a more realistic star, the results should be proportional to b 2 . 

Finally, in Ref. || we use a simple cut-off as regularization procedure, obviously different 
from both the Pauli- Villars and the dimensional scheme used, respectively, by us in Refs. 
I]!! and by Kiers and Tytgat in Ref. H; we then perform a "crude" approximation, retaining 
only the sharp effects of the neutrino potential. Our result is proportional to b 2 (Eq. fllij)). 
We might wonder why this "crude" result does not show larger powers of b. We may 
conjecture that this is due to the fact the two-body potential, being the shortest range one, 
is the one that corresponds to the "sharp" effects considered in the "crude" approximation, 
while the many-body potentials, being long-range, contribute to the corrections to the "crude 
approximation" studied in 0. 

As a conclusion we believe the result fll4]) to be quite reasonable. However, our main 
conclusion is that the energy should be even in b because of CP invariance and that the 
symmetrization given by Eq. (|9]) is required. That the resulting energy is even in b is 
also manifest in Schwinger's expansion ||, as it results from the symmetrization (Furry's 
theorem) implicitly assumed in a Feynman diagram treatment of the problem. 
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K. Kiers and M. Tytgat for very important comments about that question. This work has 
been partially supported by Spanish CICYT, project PB 95-0533-A, and by the Fundacion 
Ramon Areces. 



8 



REFERENCES 

[1] A. Abada, M.B. Gavela and O. Pene, Phys. Lett. B387 (1996) 315. 
[2] M. Kachelriess, Phys. Lett. B426 (1998) 89. 

[3] A. Abada, O. Pene and J. Rodriguez-Quintero, Phys. Lett. B423 (1998) 355. 



[4] K. Kiers and M.H.G. Tytgat, Phys. Rev. D57 (1998) 5970; |hep-ph/98074T2 . 



[5] A. Abada, O. Pene and J. Rodrfguez-Quintero, [hep-ph/980239l , Phys. Rev. D58 (1998) 
073001. 



[6] J. Arafune and Y. Mimura, |hep-ph / 9805395| , to appear in Prog. Theor. Phys. 



[7] J. Schwinger, Phys. Rev. 94 (1954) 1362. 

[8] H. Kloor, E. Fischbach, C. Talmadge and G.L. Greene, Phys. Rev. D49 (1994) 2098; 
E. Fischbach, Ann. Phys. 247 (1996) 213; B. Woodahl, M. Parry, S-J. Tu and E. Fis- 
chbach, |hep-ph/9709334| ; E. Fischbach, |hep-ph/9606250| ; E. Fischbach and B. Woodahl, 



|hep-ph/980T387 



[9] C. Itzykson and J. B. Zuber, "Quantum field theory" (McGraw-Hill, New York, 1980). 

[10] G. Feinberg and J. Sucher, Phys. Rev. 166 (1968) 1638; J.B. Hartle, Phys. Rev Dl 
(1970) 394; G. Feinberg, J. Sucher and C.-K. Au, Phys. Rep. 180 (1989) 83; 
S.D.H. Hsu and P. Sikivie, Phys. Rev. D49 (1994) 4951. 



9 



